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ABSTRACT 

We present time dependent modeling based on the accretion disk limit cycle model for a 90 d light curve of 
the short period SU UMa-type dwarf nova V344 Lyr taken by Kepler. The unprecedented precision and cadence 
(1 minute) far surpass that generally available for long term light curves. The data encompass a superoutburst, 
preceded by three normal (i.e., short) outbursts and followed by two normal outbursts. The main decay of the 
superoutburst is nearly perfectly exponential, decaying at a rate ~12 d mag -1 , while the much more rapid 
decays of the normal outbursts exhibit a faster-than-exponential shape. We show that the standard limit cycle 
model can account for the light curve, without the need for either the thermal-tidal instability or enhanced mass 
transfer. 

Subject headings: accretion, accretion disks - binaries: close - cataclysmic variables - stars: dwarf novae 


1. INTRODUCTION 

Dwarf novae constitute a subclass of the cataclysmic vari- 
ables, semi-detached interacting binaries in which a Roche- 
lobe filling secondary transfers matter to a more massive, and 
also more compact primary. The transferred material pos- 
sesses angular momentum, and therefore can only “fall” to- 
ward the primary down to a radius determined by its spe- 
cific angular momentum (Lubow & Shu 1975). This angu- 
lar momentum barrier is overcome by angular momentum 
transport within the accretion disk which carries angular mo- 
mentum outward and mass inward. Current thinking about 
the physical mechanism responsible for accretion centers on 
the magneto-rotational instability (=MRI; Balbus & Hawley 
1998) in which the shearing amplification of a weak seed 
magnetic field leads to turbulent transport. 

V344 Lyr is a short period dwarf nova, of subtype SU UMa 
(P or b = 2.06 h; Still et al. 2010) exhibiting a supercycle pat- 
tern (109.6 d; Kato et al. 2002) consisting of superoutbursts 
separated by normal outbursts. Superoutbursts are long out- 
bursts with slow, exponential decays exhibiting superhumps 
- modulations in the light curve at a period greater than the 
orbital period by a few percent. The data presented in Kato 
et al. reveal only the brightest portions of the superoutbursts; 
the fainter portions and all of die normal outbursts lie below 
the lower limit of detection. 

The Kepler data set for V344 Lyr ( Kepler ID 7659570) — 
90 d at 1 min cadence - provides an extraordinary resource 
for accretion disk modelers. The accretion disk limit cycle 
mechanism, which is currently employed to account for dwarf 
nova outbursts, is based on the storage of material during 
quiescence in a non-steady-state configuration, followed by 
a dumping of matter onto the accreting WD when a critical 
surface density is attained in the disk, producing an outburst. 
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During outburst the disk is roughly in steady state, except for 
the action of a cooling front that begins at file outer disk and 
moves inward. The thermal-viscous limit cycle model is not a 
complete model for the accretion disk, but rather primarily a 
model for the outbursts; for instance it fails to account for the 
high X-ray luminosities and column densities seen in dwarf 
novae in quiescence (Mukai, Zietsman, & Still 2009). 

Previous attempts at constraining physical parameters of 
the model have relied on AAVSO (American Association of 
Variable Star Observers) and RASNZ (Royal Astronomical 
Society of New Zealand) data for the long term light curves. 
These data are typically given as Id means, with large at- 
tendant error bars (~0.3— 0.5 mag). The Kepler data pro- 
vide an entirely new level of precision, and allow detailed 
model constraints. The viscosity parameter a (Shakura & 
Sunyaev 1973) sets the disk time scales. To first order, a co id 
controls the recurrence time for outbursts, and c* hot sets the 
duration of the outbursts. Previous efforts, for instance by 
Smak (1984=S84), utilized the “Bailey relation” relating the 
rate of the fast decay in DNe with orbital period to infer that 
Qihot — 0.1 - 0.2. This remains the only firm constraint we 
have in astronomy for a in ionized gas, and has served as 
a benchmark value for modelers calculating the efficiency of 
the MRI in accretion disks (e.g., Hirose et al 2009). 

Our goal in this work is to find model parameters that can 
account for the gross properties of the long term V344 Lyr 
light curve. There will be additional features in the light curve 
that can be compared to the detailed model light curves in or- 
der to gain insight into shortcomings of the model. We also in- 
vestigate systematic effects such as the number of grid points 
and time step size. 

In section 2 we present an overview of the physics of ac- 
cretion disks. Section 3 discusses the V344 Lyr light curve, 
section 4 presents a review of viscous decays in several well- 
studied dwarf novae, section 5 details the numerical model, 
section 6 presents the results, section 7 contains a discussion, 
and section 8 sums up. 

2. ACCRETION DISK PHYSICS 

In the accretion disk limit cycle model, gas accumulates 
in quiescence and accretes onto the central object in out- 
burst (e.g., Cannizzo 1993a, Lasota 2001 for reviews). The 
phases of quiescence and outburst are mediated by the action 
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of heating and cooling fronts that transverse the disk and bring 
about phase transitions between low and high states, consist- 
ing of neutral and ionized gas, respectively. During quies- 
cence, when the surface density £(r) at some radius within 
the disk exceeds a critical value £ ma x(r), a transition to the 
high state is initiated; during outburst, when £(r) drops be- 
low a different critical value £ m in(0, a transition to the low 
state is initiated. Low-to-high transitions can begin at any ra- 
dius, whereas high-to-low transitions begin at the outer disk 
edge. This situation comes about because in the outburst disk 
£(r) oc r -3 / 4 (roughly), and the critical surface densities 
both increase with radius. Since the disk mass accumulated 
in quiescence is bounded by SmaxM and £ m in(r), one can 
define a maximum disk mass 


J' 27rrdr£ max (v ) 


and a minimum disk mass 


Afdisk, min — J 27rrdr£ m i n (r) (2) 

which will bound the general, time dependent disk mass. 

For normal, “short” outbursts, only a few percent of the 
stored gas accretes onto the central object: the thermal time 
scale of thin disk is short compared to the viscous time scale, 
and the cooling front that is launched from the outer edge of 
the disk almost as soon as the disk enters into outburst tra- 
verses the disk and reverts it back to quiescence. For disks that 
have been “filled” to a higher level with respect to M disk) max , 
the surface density in the outer disk can significantly exceed 
the critical surface density £ m i n . In order for the cooling front 
to begin, however, the outer surface density £(r oute r ) must 
drop below £ m in(r 0 uter)* Disks in this state generate much 
longer outbursts, with slower “viscous” plateaus, because the 
entire disk must remain in its high, completely ionized state 
until enough mass has been lost onto the WD for the condi- 
tion £(r ou ter) < E m i n (r oute r) to be satisfied. The slow de- 
cay associated with a superoutburst is a direct reflection of the 
viscous time scale in the outer disk, whereas the subsequent 
faster decay reveals the thermal time scale. Although the de- 
cays of dwarf novae at a given orbital period are fairly uni- 
form, and form the basis for the Bailey relation, the rise times 
show a greater variety, reflecting the fact the the outburst can 
be triggered anywhere in the disk (Cannizzo 1998b). Inside- 
out bursts tend to produce slow-rise times, whereas outside- 
in bursts produce fast rises (Cannizzo, Wheeler, & Polidan 
1986). 

Whitehurst (1988) discovered a hydrodynamical instability 
that can occur in the outer parts of accretion disks when the 
mass ratio q = M 2 /M 1 is less than 1/4. For these extreme 
mass ratio systems, expected in CVs below the period gap if 
the secondary fills its Roche lobe, the outer part of the disk 
can lie beyond the point of 3:1 resonance with the binary pe- 
riod. The outer disk can then be excited by the inner Lindblad 
resonance (Lubow 199 lab), causing the global disk oscilla- 
tion mode that is observed as common superhumps. Before 
the onset of the superhump oscillation, there are well-known 
spiral dissipation waves in the disk that are fixed in the co- 
rotating frame. Once the superhump oscillation begins, these 
spiral arms advance 180° every superhump cycle 7 . As a given 
arm passes between the two stars, it expands outward in the 
shallower potential, but as fluid in that arm compresses back 
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Fig. 1. — Kepler light curve of V344 Lyr over 89 d showing outbursts 
every ~10 d, with one superoutburst and one failed normal burst. The flux is 
measured in e~ cadence - 1 , where the cadence (time between integrations) 
is 1 m. The full time resolution is shown. To obtain visual magnitude, a 
correspondence of my = 12 to 10 7 e~ cadence -1 was adopted. 

into the disk as the disk is maximally distorted, viscous dissi- 
pation causes a peak in the photometric light curve (Simpson 
& Wood 1998). A secondary source of the superhump signal 
is caused by the variable distance the mass stream must fall 
from the LI point to the edge of the non-axisymmetric, flex- 
ing disk. Because the outer disk is thought to expand during 
superoutburst, potentially beyond the 3:1 radius, this “tidal 
instability” model provides a natural explanation for super- 
humps. If the mass ratio is extreme enough, superhumps can 
also occur outside of superoutbursts because the disk can re- 
mained extended beyond the 3:1 radius, and therefore eccen- 
tric, even into quiescence. Osaki (1989ab) combined the ac- 
cretion disk limit cycle model with Whitehurst’s tidal instabil- 
ity model and proposed the “thermal-tidal” instability model 
(TTIM) for superoutbursts. This model posits that the pres- 
ence of an eccentric precessing disk, triggered by the migra- 
tion of gas beyond the 3:1 resonance radius, also leads to a 
greatly enhanced tidal torque acting on the disk. This ulti- 
mately contracts the outer disk and leads to a superoutburst. 
Thus in the TTIM one has the usual hysteresis relation be- 
tween surface density and temperature in the disk, and also 
a hysteresis relation between the total angular momentum of 
the disk and the tidal torque acting on the outer disk. (The 
latter hysteresis is somewhat speculative.) This model was 
examined in detail and compared to the enhanced mass trans- 
fer model (EMTM) for superoutbursts by Ichikawa, Hirose, 
& Osaki (1993) and later workers, most recently Schreiber et 
al (2004). 

3. KEPLER LIGHT CURVE OF V344 LYR 

Figure 1 shows the Kepler light curve of V344 Lyr. The ob- 
servation window contains outbursts showing the sequencing 
SSSLSS, where S= short outburst and L=long outburst, i.e., 
a superoutburst. Based on photographic observations made 
by Hoflmeister (1966), Kato (1993) estimated that the cycle 
length for normal outbursts was 16=b3 d. An estimate of this 
cycle length is difficult because the normal outbursts of V344 
Lyr are short-lived and reach only to near my = 16. Using 
our five complete normal outbursts, we find that prior to the 



7 See, for example, www.astro.fit.edu/wood/visualizations.html. 
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TABLE 1 

Duration and Amplitude of Normal 
Outbursts 


Outburst 

Duration (day) 

Amplitude (mag) 

1 

5.5 

4.0 

2 

5.0 

3.6 

3 

5.3 

3.6 

4 

3.3 

2.4 

5 

4.0 

2.6 


superoutburst, the cycle length is 9.6 d, while after the super- 
outburst it shortens to 7.2 d. The normal outburst amplitude 
and duration also changes from before (normal outbursts 1 , 2, 
3) to after (normal outbursts 4, 5) the superoutburst (Table 1). 
We also note that the level of quiescence is ~0.5 mag higher 
after the superoutburst, indicating a larger contribution from 
the disk emission. Since the disk is a bit smaller and has less 
mass at this point (see Section 6), we believe that this brighter 
disk is due to it being hotter overall directly after the super- 
outburst. 

The superoutburst lasted 17.6 d and reached a maximum 
amplitude of 4.2 mag. The decline of the superoutburst is 
close to linear (plotted as mag versus time, or exponential as 
flux versus time) from near day 57 until day 68, after which 
it falls rapidly (~2 d) back to quiescence. During this part 
the decline rate is 0.083 mag d 1 , in good agreement with the 
value of 0.094 mag d -1 found by Kato (1993). The amount of 
time spent at my < 16 is ~ 14 d, and the duration of the entire 
plateau is ~12 d. The decay rate is ~12 d mag -1 , therefore 
the plateau portion of the decay encompasses ~1 mag — a 
dynamic range of ~2.5, or slightly less than one e — folding 
(~2.7x). There is also a ~0.5 mag shoulder on the superout- 
burst rise lasting rsj 1 

It is noteworthy that the fast decays of the normal out- 
bursts and superoutburst have a faster-than-exponential shape 
(i.e., concave downward when plotted as magnitude versus 
time). S84 utilized the Bailey relation between the fast rate 
of decay in dwarf novae and their orbital period to constrain 
a hot — 0.1-0.2. The large errors associated with the AAVSO 
data were consistent with exponential decay. Cannizzo (1994) 
found that to account for the supposedly exponential decays, 
which correspond to the time during which a cooling front 
traverses the disk, a must vary weakly with radius (~oc r 0 3 ). 
As we show later, the faster-than-exponential decays captured 
by Kepler in V344 Lyr negate this result and indicate that a 
constant a scaling with radius is consistent with the data. 

For any dynamical system in which a variety of timescales 
contribute to a physical observable, in this case the flux from 
the accretion disk, the controlling time scale will be the slow- 
est one. Since all the disk timescales — dynamical, thermal, 
and viscous — increase with radius, the time scale associated 
with the viscous plateau provides direct information about the 
viscous time scale t v isc at the outer edge of the disk. 

Starting with a general theoretical expression from Frank, 
King, & Raine (2002) relating the locally defined vis- 
cous time to a , r, and accretion rate M, King & Pringle 
(2009) estimate a local peak accretion rate during outburst 
as ~2M disk; hot Aviso where M diskj hot is the mass of the hot 
disk which they estimate as M diskj mill times a factor of order 
unity, and derive a viscous time for the outbursting disk (see 
their eqn. [2]) 

*visc = 14 d oT?' 8 r? o 48 m? 15 (3) 


where a_i = a/0.1, ri 0 = r/(10 10 cm), and mi = 
Mi/M©. (The a in this context is the hot state value.) Within 
the viscous plateau portion of a superoutburst, the instiga- 
tion and propagation of a cooling front at the outer edge is 
thwarted by virtue of excessive mass in the disk. The viscous 
time does not represent the duration of the outburst itself, but 
rather the time scale (approximately in d mag -1 ) associated 
with the viscous plateau; it represents an e— folding time for 
the decrease of disk mass. 

4. VISCOUS PLATEAUS IN DWARF NOVAE 

The Bailey relation between the fast rate of decay of dwarf 
nova outbursts and orbital period corresponds to the time in- 
terval during which the cooling front traverses the disk. Some 
long dwarf nova outbursts have enough dynamic range in my 
that one can also use the slow decay portion, the viscous 
plateau, to infer an a value. These inferred values agree with 
those determined from the Bailey relation (Cannizzo 2001b, 
Cannizzo et al. 2002; King et al. 2007). 

The decay rates associated with several dwarf novae ex- 
hibiting viscous plateaus have been fairly well established, 
and enable one to make a reasonable estimate of the decay 
rate expected in a system with the orbital period of V344 Lyr. 
As noted previously, for V344 Lyr one expects t visc ~ 14 
d mag -1 . The Kepler observation gives a decay rate within 
the superoutburst ~12 d mag -1 . We associate this with the 
viscous time, which is basically an e-folding time for a per- 
turbation to the £(r) profile to be damped out. For a disk with 
a £(r) profile approximating steady state, i.e., £(r) oc r -3 / 4 , 
and for which the mass flow within the disk and onto the 
WD greatly exceeds the rate of mass addition feeding into the 
outer disk from the secondary star, the viscous time represents 
roughly a time for the mass of the disk to decrease by a factor 
e. 

Thus for a given viscous plateau-type outburst, the ratio 
^plateau / ^visc represents the number of e— foldings by which 
the mass of the disk decreases due to accretion onto the central 
object. It is instructive to look at examples of viscous decays 
in several well-studied dwarf novae: 

SS Cyg: Long outbursts in SS Cyg (P orb = 6.6 hr) have a 
duration of ~10 d (Cannizzo & Mattei 1992). and the viscous 
time in the outer disk £ visc » 40 d. Thus t p i a teau < t v isc» and 
only a small fraction of the disk mass ~0.2 is accreted onto 
the WD (Cannizzo 1993b=C93b). 

U Gem: In the 1985 long outburst of U Gem (P orb = 4.25 
hr), the decay rate during the viscous plateau was ~26 d 
mag -1 over the ~35 d of the burst (Cannizzo et al. 2002). 
Since t v i sc « plateau, the disk mass decreased to about 0.3 of 
its initial value. 

V344 Lyr: In the Aug 2009 superoutburst of V344 Lyr 
(P orb = 2.06 hr) shown in Fig. 1 the locally defined decay 
rate during the plateau remained constant at ~12 d mag -1 
over most of the ~ 14 d of the burst, thus as with U Gem the 
disk mass remaining after the superoutburst would be about 
0.3 of the initial amount. 

WZ Sge: In the 2001 superoutburst of WZ Sge (P orb = 81 
min), the locally defined decay rate during the plateau in- 
creased from ~ 4 to ~ 12 d mag -1 over the ~20 d of the burst 
(Cannizzo 2001b). The fact that t p i a teau ~ 3t v i sc means that 
only a few percent of the initial disk mass remained at the end 
of the superoutburst. (An alternative considered in the Discus- 
sion is that irradiation-induced enhanced mass overflow from 
the secondary star prolonged the superoutburst.) 

In summary, one can readily find examples of known sys- 



4 


CANNIZZO ET AL 


terns with viscous type outbursts in which the duration of the 
outburst itself is either greater than, less than, or about the 
same as the viscous time scale evaluated at the outer edge of 
the accretion disk. The actual course taken by a given system 
depends on the degree of overfilling of the disk in quiescence 
with respect to M disk5min . 


Some early studies adopted C = D = 0 (M83), C = 0 (L85), 
or a form of C not expressible as the divergence of a flux 
(M83, M86, C93b; see H98 for a critical discussion). 

The local flow velocity v r , which can vary enormously from 
inflow to outflow over the radial extent of the disk, depending 
on whether or not transition fronts are present, is given by 


5. NUMERICAL MODELING 

Our numerical model is discussed extensively in C93b and 
subsequent papers. The basic strategy is the time depen- 
dent solution of two coupled differential equations, one giving 
the £(r) evolution, and one the evolution of disk midplane 
temperature T(r). Thus the disk is entirely unconstrained 
as regards deviations from steady state and thermal equilib- 
rium. Many previous workers have presented time depen- 
dent calculations by solving these equations, or similar ver- 
sions (e.g., S84; Lin, Papaloizou, & Faulkner 1985=L85; Mi- 
neshige & Osaki 1983, 1985=M85, Mineshige 1986=M86, 
1987=M87; Pringle, Verbunt & Wade 1986=P86; Cannizzo 
et al. 1986=C86, C93b, 1994, 1998a, 2001a; Angelini & 
Verbunt 1989; Ichikawa & Osaki 1992; Ludwig & Meyer 
1998; Hameury et al. 1998=H98, Menou, Hameury, & Stehle 
1999, Buat-Menard, Hameury, & Lasota 2001, Buat-Menard, 
& Hameury 2002, Schreiber, Hameury, & Lasota 2004= S04). 

Improvements in our code made since the original version 
include the provision for a variable outer disk radius, self- 
irradiation, and evaporation with various possible laws. The 
diffusion equation governing the evolution of surface density 
is given by (Lightman 1974, Pringle 1981) 



where £(r,t) is the surface density and v is the kinematic 
viscosity coefficient, 


2 a 1ZT 

3 D fi 


(5) 


We follow the technique of Bath & Pringle (1981) in discretiz- 
ing the £(r, t) evolution equation. 

The thermal energy equation governing the evolution of the 
midplane disk temperature T(r , t) is given by 


dT 2 (A-B + C + D) IZTld, x dT ^ 

at ^ — 

where the viscous heating 


the radiative cooling 


A = -iyQ 2 E, 


B = oTt, 


r= 31_9 / 

2 r dr \ 


^ r %) 


and 


„ h d f 4 acT 3 dT' 
D = — — r- 


f 4 acT 3 dT\ 
r dr \ 3krp dr J 


(7) 

( 8 ) 
(9) 

( 10 ) 


Of the four terms appearing in the numerator of the first 
term of thermal energy equation, the first two A and B are 
fairly standard. The term C represents the radial heating flux 
due to turbulent transport, and D is that due to radiative trans- 
port. Previous workers have used varying forms of C and D. 


O') 

where upwind differencing is implemented over the radial 
grid as follows: 


s a = i/i_i£i_iri_i 1/2 

s b = I 'i5W /2 

Sc = ^i+l^i+lTi+1 ^ (12) 


8 Sf) s a 

EiJ-j 1 / 2 Tj - V{— l 


3 s c Sjj 

Sin 1 / 2 r i+ i - n 


(13) 


if 

if 


— >0 and q- > 
9+ 

— > 0 and q- < 
Q+ 


°) 

°) 


(«r)i =max(g_,<jf+) 


K)i =min(<jf_,<ji + )(14) 


if I — < 0 and j^4 > 1^ 


9+ 


\q+\ 


it| ^<0a»d t}<l' 
9+ 19+1 j 


K)i =9- 
(v r )i=q+. (15) 


The scalings for the local maximum and minimum in the 
£ — T relation (from MM81) 


Smax = 444 gem 2 r\^m 1 


0 . 37 - 0.7 


OL, 


C,~2 5 


(16) 


where a C5 -2 = a C oid/0-C)l> and 

Smin = 44.4 g cm -2 °- 37 a^i, (17) 


where a^-i = Qihot/0.1. The disk midplane temperatures 
associated with these extrema are 


T(S max ) = 5275Ka-°_ 2 3 , (18) 

and 

T( E min ) = 35900 K rX^m^a^ 6 . (19) 

The equilibrium temperature scalings for the stable branches 
are given in C93b (taken from Cannizzo & Wheeler 1984). 
During transitions when the local disk temperature lies be- 
tween the values associated with the local maximum and min- 
imum in surface density, T[£ max (r)] < T(r) < T[£ mill (r)], 
we use a logarithmic interpolation to obtain the local a value, 

logi 0 Q!(r) = log 10 a cold + /, (20) 

where 

/ = ^-^^(logiO^t - logiottcold)- (21) 
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In tests we also consider a logarithmic interpolation factor / 
introduced by H98, 


f = 




(l°glO a hot l°g 10 ^cold ) 5 


( 22 ) 


as well as a linear interpolation 


ot(r) = a co id + 


T-T{Y max ) 


(oihot ^cold)- (23) 


In calculating the effects of the tidal torque from the sec- 
ondary we follow Smak (1984) and later workers who used 
a formalism developed by Papaloizou & Pringle (1977) in 
which the tidal torque varies as the fifth power of radius. This 
leads to a depletion in disk material 


dY vY /r\ 5 

dt ClUJ 2i vj s Va/ 


(24) 


where u is the orbital angular frequency 27r/P 0 rb> the or- 
bital separation is a, and the specific angular momentum 
j s = (GMir) 1 / 2 . We do not reset the value of ci during the 
course of a run as in the TTIM. Ichikawa & Osaki (1994) pre- 
sented detailed calculations of the strength of the tidal torque 
versus distance at large radii within the disk and found that 
the power law (oc r 5 ) determined by Papaloizou & Pringle is 
only valid at one radial location in the outer disk; further out 
the torque relation steepens considerably, approaching infi- 
nite slope at the last non-intersecting orbit. Therefore in some 
sense the more naive “brick wall” condition for the outer edge 
used by C93b and other workers might be more descriptive. 
Regardless of which prescription is used, if the parameters 
associated with the outer edge are not manipulated during 
the run, the detailed treatment is not important to the over- 
all model. Using a different end point treatment would result 
in slightly different optimal parameters. 

Given a disk temperature T and surface density E, the disk 
semithickness h and density p are determined from vertical 
hydrostatic equilibrium to yield 


u P r + y/Pp + C2C3 

^ ~ ~ 5 
C2 

(25) 

where 

C2 = EJ2 2 , 

(26) 

, y nT 
c 3 — Y , 

(27) 

the radiation pressure 


Pr = \aT\ 

(28) 

and the density 

E 

P ~ 2ft' 

(29) 


To determine the computational time step, we take 


At = ft min 


min; 




, min; 




(30) 

where min; refers to the minimum over all grid points i, and 
ft is a small number between 1/80 and 1/640. The quantities 
dlnY i/dt and 9 In Ti/dt are evaluated using the evolution 
equations. 
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Fig. 2. — The effect of using a logarithmic interpolation (eqn. [21]) versus 
a linear one (eqn. [23]) for a. For these runs N = 800 and ft = 1/160. 
Shown are the light curves for the logarithmic case ( top panel) and the linear 
case {second panel), as well as the accompanying disk masses {third and 
fourth panels ). The disk masses are in units of 10 22 g. 
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Fig. 3. — The effect of using the H98 logarithmic factor / (eqn. [22]). 
For these runs N = 800 and ft = 1/160. Shown are the light curves 
for To = 25000 K {top panel), 18750 K {second panel), and 12500 K 
{third panel), as well as die accompanying disk masses (fourth through sixth 
panels ). The disk masses are in units of 10 22 g. 


6. RESULTS 

To calculate the long-term light curves from our numerical 
trials, we assume local Planckian flux distributions for each 
annulus, and in each time step integrate over the radial profiles 
of effective temperature. We compute the absolute Johnson 
V -band magnitude, assuming a face-on disk. A reasonable 
inclination i ~ 45° would dim the disk by ~0.4 mag. After 
some experimentation, we found a set of parameters that re- 
produces approximately the overall pattern of supercycle and 
short outbursts. We take a central O.6M 0 WD accretor, an in- 
ner disk radius r in = 2 x 10 9 cm, and a maximum outer disk 
radius r ou t = 10 10 cm, a mass feeding rate into the outer disk 
M t = 2 x 10 _1 °M o yr -1 , and a values a co \d = 0.0025 
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Fig. 4. — Shown are the different a co ia — cxhot interpolations. The 
blue lines show the scalings for T(£ max ) and T(E min ). The black curves 
indicate the linear (eqn. [23]) and logarithmic (eqn. [21]) interpolations. The 
green curves show the H98 scaling (eqn. [22]), where the three curves (left to 
right) represent To = 25000 K, 18750 K, and 12500 K. The three red curves 
show a measure of the partial ionization; they are the relative contribution 
of electron pressure to total pressure as shown in Fig. 16 of C93b, for p ( g 
cm -3 )= 10 -8 , 10 -7 , and 10“ 6 (left to right). The values given in C93b 
have been multiplied by 2o:hot t0 take into account the fact that (1) P e /P —3 ► 
0.5 in the limit of complete ionization, and (2) in the high temperature limit 
OL tthot = 0.1. 

and ahot = 0.1. We use the a interpolation given by eqn. 
(21). As regards the inner disk radius, the value for r inner is 
larger than the WD itself; if we take ri nner ~ rw d, we tend 
to get far fewer short outbursts between two successive super- 
outbursts. This hint of a large inner radius may have some 
observational support in the SU UMa systems (e.g., Howell 
et al. 1999). The ratio a C oid/^hot = 40 is about 10 times 
greater than commonly used in dwarf nova calculations. The 
mass transfer rate from the secondary M T is twice the nom- 
inal value 10 _10 M© yr -1 found from evolutionary calcula- 
tions of cataclysmic variables below the period gap where the 
evolution is driven predominantly by gravitational radiation 
(e.g., Kolb, King, & Ritter 1998; Howell, Nelson, & Rappa- 
port 2001). When one takes into account the ~20% deple- 
tion of the disk mass accompanying each of the short (nor- 
mal) outbursts, however, the effective rate of accumulation of 
mass in the disk between successive superoutbursts decreases 
from 2 x 1O -1O M 0 yr -1 to ~0.3xl0 _1 °M o yr -1 . In our 
models, ~40-60% of the disk mass is accreted during a su- 
peroutburst. 

In using a complicated numerical model, it is important to 
carry out testing to gain an understanding of systematic ef- 
fects: 


6.1. Interpolation between a C oid and ahot 

Figure 2 shows the effect of using a log versus a linear in- 
terpolation for the a value between a co id and a^ot . The linear 
interpolation leads to smaller amplitude outbursts with little or 
no quiescent intervals. The decays are also much slower, and 
concave upwards when plotted as mag versus time (i.e., with 
a functional form that is slower-than-exponential). Figure 3 
shows the effect of using the H98 interpolation. H98 adopted 
the constant T 0 = 25000 K in their scaling; we also look at 
smaller values 18750 K and 12500 K. The smaller T 0 values 
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Fig. 5. — Long-term model light curves for N = 100, 200, 400, 800, and 
1600. For these runs ft = 1/ 640. The light curves have been aligned so that 
a superoutburst begins at t = 25 d. 
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Fig. 6. — The disk mass in units of 10 22 g accompanying the light curves 
shown in Fig. 2. 

produce superoutbursts which have the recognizable plateau 
and subsequent fast decay, but for the T 0 = 12500 K run the 
fast decays are slower-than-exponential. 

Figure 4 depicts the different interpolations as well as pa- 
rameterizations of the fractional ionization £. The MRI is 
thought to mediate the strength of a through turbulent trans- 
port, and the effect of shearing on the weak magnetic field 
embedded in the gas is more efficient as £ increases. The H98 
scaling for T 0 = 12500 K does a reasonable job in approxi- 
mating the p = 10 -6 g cm -3 curve for £, and may be more 
physical than the other scalings. 

6.2. Number of Grid Points N 

Early time dependent accretion disk limit cycle studies were 
concerned mainly with showing the ability of the disk, under 
the conditions of an imposed limit cycle in each annulus, to 
undergo collective, global oscillations in which a major part 
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Fig. 7. — Model light curves showing the detail of a superoutburst from 
each panel of Fig. 5 for N = 100 {dotted), 200 {short dashed ), 400 {long 
dashed ), 800 {short dash-dotted ), and 1600 {long dash-dotted). The green 
curve shows a superoutburst from the T 0 = 12500 H98 scaling light curve 
of Fig. 3. For comparison the averaged Kepler data for V344 Lyr is shown 
(blue), where the original At = 1 min cadence has been block-averaged to 
At = 4 hr. The two slanted red lines indicate a decay slope of 12 d mag -1 . 
The second panel shows the disk mass (in 10 22 g) as in Fig. 3. 

of the disk participated. The total numbers of grid points N 
was small and the comparison with observational data mini- 
mal (e.g., N = 20 - 25: S84; N = 35: L85; N = 22 - 44: 
M85; N = 17: M86; N = 40: P86; N = 44: C86). One can 
see obvious irregularities in many of the early light curves 
which are characteristic of numerical instability. C93b un- 
dertook the first detailed investigation regarding the number 
of grid points, and found that at least several hundred were 
required to obtain reliable light curves. H98 utilized both uni- 
form and adaptive mesh codes to study the properties of the 
transition fronts, and found a requirement of N > 800 for a 
uniform-equivalent mesh to resolve the fronts adequately. 

Figure 5 shows long term light curves for models with 
N = 100 to 1600 grid points. The rise times we calculate, 
both for the normal and superoutbursts, are too fast. In ad- 
dition, our outbursts are too faint by ~0.5 mag. A nominal 
inclination i 45° would make our disks dimmer by ~0.4 
mag and increase the discrepancy to ~1 mag. The shoulders 
on the rise to superoutburst in the models have more structure 
than the shoulder in the observed superoutburst. The shoul- 
ders are caused by the strong spreading of material out of the 
annulus in which the instability starts. This results in a brief 
reversal in the increase of the local disk temperature. The rate 
of brightening in V is set by a combination of the widening 
radial extent of the unstable annulus, and the increase or de- 
crease of temperature within the annulus. The subsequent rise 
after the shoulder is due to the full radial extent of the disk out 
to the outer edge making a transition to the hot state. In con- 
trast, for the short (normal) outbursts, not enough disk mass 
has accumulated at larger radii by the time of the outburst to 
support the heating transition front propagating to large radii, 
therefore the heating front stalls at some r < r ou ter and is 
reflected as a cooling front, leading to pointy maxima rather 
than viscous plateaus. 


For smaller N one sees spuriously long and slowly decay- 
ing superoutbursts, and only 3 short outbursts between two 
successive superoutbursts. For N > 800 the pattern asymp- 
totes to LSSSSSSSLSSSSSSS... The decrease in superout- 
burst duration from N = 100 to N = 200 is consistent with 
the trend noted by C93b. Figure 6 indicates the variations in 
accretion disk mass accompanying each panel in Fig. 5. Fig- 
ure 7 shows a detail of one superoutburst from each panel in 
Fig. 5. For comparison we also plot the V344 Lyr data, where 
D = 620 pc was assumed in converting to absolute magni- 
tude (found by Ak et al. 2008 using a period-luminosity rela- 
tion). There is a slight increase in superoutburst duration from 
N = 200 to TV = 800 due to the fact that more mass accu- 
mulates as the supercycle lengthens. In addition, the shoulder 
in the light curve also becomes longer. For the highest N run 
the duration of the superoutburst is still shorter than that seen 
in V344 Lyr, while the shoulder is too long. The superout- 
burst calculated using the H98 scaling with T 0 = 12500 K 
better reproduces the observed duration of the V344 Lyr su- 
peroutburst, and also has a shorter shoulder with less structure 
than the other calculations. However, it suffers from a rise 
time that is too short, and the fast decay is now slower-than- 
exponential. Also, the viscous decay rate, which has been 
optimized to the observed value using the eqn. (21) scaling, 
is now too fast. 

6.3. Thermal Energy Equation 

Considering the many forms of the thermal energy equa- 
tion used by previous workers, and the potential limitations of 
employing a one-dimensional hydrodynamical formalism, we 
now look at the effects of switching off various terms in the 
thermal energy equation in an attempt to ascertain the relative 
importance of each term. Figure 8 shows the effect of setting 
equal to zero: (i) (7, (ii) D, (iii) C and D, (iv) the first ad- 
vective term in eqn. (6), and (v) the second advective term in 
eqn. (6). In the C = 0 run one sees a rounded superoutburst 
and only three short outbursts between two superoutbursts. 
The (radiative) D term is less influential, as the superoutburst 
profile is relatively unchanged from the first panel. Setting 
C = D = 0 produces a single dip in the rise to superout- 
burst, and a longer viscous plateau (i.e., exponential decay). 
The first advective term appears to be of comparably minimal 
importance as the D term; i.e., the omission of either has a 
negligible effect. The second advective term is the most influ- 
ential of all terms examined: without it one has much longer 
superoutbursts - about 70% of the accumulated disk mass is 
accreted, which is about the same as in V344 Lyr. However, 
the viscous plateau still comprises less than half of the super- 
outburst, rather than the > 90% seen in V344 Lyr. 

6.4. Detuning of Optimal Parameters 

There is a complicated interplay between the different pa- 
rameters that enter into the overall results. In an effort to un- 
derstand the effect each parameter has, we look at the effect of 
introducing a slow change into a co id> ol h Q t» and Mr. Figure 
9 shows the resultant variations in the light curves and disk 
masses. The primary effect of increasing these parameters is 
to shorten the recurrence times. Figure 10 shows the effect 
of varying the inner disk radius. For smaller inner disk radii, 
one has a larger dynamic range in and as a result 

there tends to be a more complete emptying of the disk dur- 
ing outbursts, and therefore less of a tendency to have a series 
of short outbursts leading up to a superoutburst. The net in- 
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Fig. 8. — The effect of turning off different terms in the energy equa- 
tion. For these runs TV = 800 and ft = 1/160. Shown are light curves 
using: the full equation (top panel), die turbulent transport C = 0 (second 
panel), the radiative transport D = 0 (third panel), both (7 = 0 and D = 0 
(fourth panel), the first advective term set to zero (fifth panel), and the second 
advective term set to zero (sixth panel). 

teraction among all the model parameters is such that the ob- 
served sequencing for V344 Lyr is roughly reproduced with 
Tinner — 2 x 10 9 cm. Figure 1 1 shows the effect of varying the 
outer disk radius. For larger r oute r values not only does the 
number of short outbursts increase, but also the shape of the 
superoutburst distorts significantly from the viscous plateau 
shape. 

6.5. Relevance of Results to the V344 Lyr light curve 

No one light curve we have shown reproduces all features 
of the observed light curve shown in Fig. 1. For instance, 
there is an atypically long quiescent period at the beginning 
of the observed light curve. Figure 9 in which we vary a co \d> 
a hot, and Mr shows that if any one of these parameters were 
anomalously small, that could delay an outburst. Also, the 
pre-superoutburst normal outbursts are brighter than the post- 
superoutburst normal outbursts, and the pre-superoutburst 
quiescent level is fainter than the post-superoutburst quies- 
cent level. As noted earlier, the disk instability model is ill- 
equipped to make any statement about the quiescent state, but 
Figure 9 again shows several possibilities. For instance, if M T 
were elevated for a while after the superoutburst, that could 
lead to crowding of the short outbursts, and an increase in the 
quiescent level. There may be other causal relations among 
G! cold * ahot> and M t , which are important. 

The primary interpolation scheme we considered, given 
in eqn. (21), has clear problems in making good superout- 
bursts. The model superoutbursts have shoulders which com- 
prise about a third to half of the total viscous plateau, rather 
than the ~10% observed. However, for the H98 scaling with 
T 0 = 12500, the shoulders are appropriately short, ~2 d, and 
the viscous plateaus much longer. The rate of decay of the 
viscous plateau is too fast, but that is because the input pa- 
rameters have been optimized for the scaling given in eqn. 
(21). 




t(d) 

Fig. 9. — The effect of introducing a slow change into a given parameter, 
increasing the value linearly from half to twice its optimal value. For these 
runs TV = 800 and ft = 1/160. The six panels show the resultant light 
curves for a co id ( top panel), a hot (second panel), and M T (third panel). 
Also shown are the accompanying variations in disk mass, again in units of 
10 22 g (fourth through sixth panels). 
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Fig. 10. — The effect of the inner disk radius. For these runs TV = 800 and 
ft = 1/160. The four panels show the resultant light curves for r inn er = 
2 x 10 9 cm (top panel), 1.5 x 10 9 cm (second panel), 1 x 10 9 cm (third 
panel), and 0.5 x 10 9 cm (fourth panel). 


7. DISCUSSION 

In contrast to the thermal-tidal instability of Osaki 
(Ichikawa & Osaki 1992, 1994; Ichikawa, Hirose, & Osaki 
1993) in which the tidal torque is artificially enhanced by a 
factor of 20 when the outer disk expands beyond the point of 
3:1 resonance with the orbital period, our detailed modeling 
shows that the regular and superoutbursts can produced with 
a constant set of parameters. In retrospect, one may legiti- 
mately question the motivation for the increased tidal torque 
accompanying the expansion of the disk, considering that the 
ultimate driving force behind the tidal torque is the gravita- 
tional field of the secondary, which remains unchanged. We 
argue that the superoutbursts seen in the SU UMa systems 
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Fig. 11. — The effect of the outer disk radius. For these runs N = 800 
and ft = 1/160. The six panels show the resultant light curves for runner = 
0.8 x 10 10 cm (top panel), 1.0 x 10 10 cm (second panel), 1.25 x 10 10 cm 
(third panel), 1.5 X 10 10 cm (fourth panel), 1.75 X 10 10 cm (fifth panel), 
and 2.0 x 10 10 cm (sixth panel). 

may be viewed as an extension of the long outbursts seen in 
DNe above the period gap, such as U Gem and SS Cyg (Can- 
nizzo & Mattei 1992, 1998; C93b). Continued monitoring 
of V344 Lyr should reveal the stability of the long term se- 
quencing, i.e., variations in the number of shorts sandwiched 
between two longs. 

Our testing has delineated a range of allowable parameters: 

(1) The logarithmic interpolation for a produces well- 
separated outbursts, with concave downward decays, as ob- 
served. The linear interpolation produces outbursts with min- 
imal quiescence and concave upward decays. Given the phys- 
ical underpinnings of the MRI, the logarithmic interpolation 
is better motivated since the MRI strength should depend on 
the partial ionization fraction, which is a strongly increasing 
function of temperature. 

In terms of the superoutburst, the H98 scaling with T 0 = 
12500 K produces the best overall results. The shoulder fea- 
ture only persists ~ 2 d, as observed, and the viscous plateau 
is longer. About 70% of the stored mass is accreted. The 
main deficiencies are that the fast decay (i.e., that following 
the viscous plateau) is slower-than-exponential, and the rise 
times are too fast (as with the other calculations). 

(2) The minimum number of grid points necessary is > 10 3 
(in accord with C93b and H98). Even for high AT, the se- 
quencing is not perfectly stable. 

(3) Of the terms entering into the thermal energy equation, 
the radial viscous transport term C and the second advective 
term, v r (dT/dr) appear to be the most important. The radial 
radiative flux term D and the first advective term, involving 
d(rv r )/dr , are of lesser importance. 

(4) The values for a co id = 0.0025, ahot = 0.1, and M T = 
2 x 10 10 M o yr -1 are mandated by the observed spacing and 
duration of short outbursts and the superoutburst. Increasing 


their values decreases the outburst time scales. 

(5) The values ri nner and r oute r strongly mediate the num- 
ber of short outbursts between two superoutbursts. For N = 
800 and f t = 1/160 we have ~5 shorts between two longs for 
t inner = 2 x 10 9 cm, or ~3# W d for M W d = O.6M 0 . Ob- 
servations by Howell et al. (1999) seem to favor an evacuated 
inner disk. Several possibilities are discussed by these work- 
ers, including the possibility of weak polars in the SU UMa 
systems. The outbursts in our models are triggered at inter- 
mediate disk radii, therefore the outer disk radius has mini- 
mal effect on the outburst recurrence times and durations. In- 
creasing r ou ter from 1 x 10 10 cm to 2 x 10 10 cm increases 
the number of shorts between to successive longs, and also 
distorts the profile of the superoutburst. 

Menou (2000) attempted a completely self-consistent cal- 
culation in which the local value of the fractional ionization 
parameter £ (determined from density and temperature) set 
the level of the MRI instability. The model started with a 
transition to outburst state, and then continued into outburst 
and back to quiescence. When the disk started to cool, a run- 
away transition occurred such that the disk cooled to abso- 
lute zero, effectively terminating the simulation. This came 
about basically because £ depends exponentially on temper- 
ature, through the Saha equation, whereas the level of turbu- 
lent transport, determined from the MRI, depends linearly. It 
has been recognized that, due to the very low £ values near 
disk midplane in quiescence in the standard models, the stan- 
dard S- curve relation between E and T is probably a conve- 
nient approximation for a much more complicated situation. 
In other words, there must be some additional angular mo- 
mentum transport mechanism operating in quiescence. Alter- 
natively, the quiescent disk may be essentially inviscid, and it 
may be the onset of finite and large partial ionization near the 
midplane, due to the accumulated effect of irradiation on the 
quiescent disk that is in fact the trigger for dwarf nova out- 
bursts (Gammie & Menou 1998). This is sometimes referred 
to as the “dead zone” model (Gammie 1996) to encapsulate 
the fact that the quiescent disk may consist of an “active”, irra- 
diated layer with some radial mass flow that sandwiches a rel- 
atively inert inner core. In view of these considerations, and 
the current lack of predictive power in these type of models, 
we view our adoption of a small Q! co id for an SU UMa system 
as a necessary expedient which may have little physical sig- 
nificance. It is more noteworthy that a!hot> which is dictated 
by the outburst properties, is in line with that determined from 
dwarf novae in general (King, Pringle, & Livio 2007). This is 
consistent with the notion that the MRI becomes fully active 
in outburst disks where £ » 1 and reaches a well-defined sat- 
uration level, equivalent to ~10% of the gas pressure, in the 
a— disk parlance (Balbus & Papaloizou 1999). 

With c^hot = 0.1, we match the rate of decay of the viscous 
portion of the superoutburst, and also find a ~0.5 mag shoul- 
der on the rise to superoutburst. It is not clear why previous 
detailed time dependent modelers did not find this shoulder 
(e.g., S04). The main deficiency in the models we present 
is our failure to obtain a sufficiently long viscous plateau on 
the superoutburst. The observed plateau has a decay rate of 
~12 d mag -1 and lasts for ~12 d, which corresponds to ~1 
e— folding in the decay of the disk mass, meaning about 60% 
of the disk mass should be accreted. Our viscous plateau lasts 
~6 — 8 d, with proportionally less mass accreted. The values 
for a co id, <*hot, and M T have already been optimized to get 
the supercycle and normal cycle periods of ~ 1 00 d and ~ 1 0 d, 
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respectively. One suspects that the scalings for £ ma x(acoid) 
and SminKot) are too simplistic and fail to take into ac- 
count the actual disk physics, particularly that of the quies- 
cent state, which is presently unknown. One could envision 
introducing an ad hoc multiplicative parameter for E max , for 
instance to artificially enhance it, but then one would have to 
re-adjust all the other parameters, in particular increasing M T 
which is already uncomfortably large. In view of these con- 
siderations we consider our experiment moderately success- 
ful in that we are able to reproduce the duration of the viscous 
plateau within a factor of two in an over-constrained system 
for which a co id and a hot have already been determined by the 
outburst properties. 

What about having enhanced mass transfer from the sec- 
ondary prolong the duration of the superoutburst? Smak 
(2004) presented observational evidence for enhanced hot 
spot brightness near superoutburst maxima which he inter- 
preted as being due to enhanced mass transfer brought about 
by the strong irradiation of the secondary. Osaki & Meyer 
(2003) presented a theoretical study of the effects of irradia- 
tion on the photosphere of the secondary star in SU UMa sys- 
tems and concluded that mass transfer could not be enhanced 
by a significant amount. S04 presented detailed calculations 
of both the TTIM and the EMTM, and in the end concluded 
that the EMTM was favored. Their arguments are rather gen- 
eral, however, and they state "... we have not proven the 
EMTM to be correct nor the TTIM not to work.” Both models 
have a number of free parameters in addition to those already 
present in the standard limit cycle model. For the TTIM one 
has the time delay between the attainment of the 3: 1 radius by 
the outer disk edge and the onset of enhanced tidal torque, in 
addition to the functional form and amplitude associated with 
the increase. Similarly the EMTM has as free parameters the 
degree of augmentation of M T during superoutburst, as well 
as its functional form with time. 

Our calculations, which are based solely on the standard 
limit cycle model, give alternating series of many short out- 
bursts followed by a superoutburst, but the superoutburst du- 
rations are too short. This would appears to leave open the 
possibility for enhanced mass transfer, but we are unable to 
give it a strong endorsement for several reasons. The primary 
reason is that the a interpolation is uncertain, given that we do 
not currently have an ab initio theory for a based on the MRI, 
and this weak point in the theory leads to a larger uncertainty 
in the light curves than the other systematic effects. The MRI 
is a powerful instability for which the most rapidly growing 
wavenumbers in a thin disk have growth rates of ~ 0.75ft, 
a factor of ~10 6 in < 3 rotation periods (Balbus & Hawley 
1991). One rotation period in the outer disk is 27rft _1 ~ 700 s 
r% 2 , which is much faster than any time scale associated with 
the outbursts. Therefore “a” should adjust itself rapidly to a 
value in accord with the local physical conditions (i.e., Q in 
the disk at any given time. The H98 scaling with T 0 = 12500 
may therefore represent the best physically motivated of the 
interpolation scalings we investigate, and it also does the best 
job at reproducing the main features of the superoutburst - its 
long duration, and the short (~2 d) shoulder feature. 

In addition, there may be a problem with the viscous de- 
cay in the EMTM: The V344 Lyr superoutburst maintained 
a closely exponential decay over ~1 mag, whereas our mod- 
els as well as others reveal a definite slower-than-exponential 
trend (i.e., concave upward when plotted as magnitude versus 


time) in viscous plateaus covering that much dynamic range. 
Although the viscous decays can be lengthened in the EMTM, 
they tend to have noticeable deviations from exponentiality, 
in contrast to the V344 superoutburst. The fact that the short 
orbital period system WZ Sge did show a deviation from ex- 
ponentiality in its 2001 outburst, which also spanned a greater 
dynamic range than in the V344 Lyr superoutburst (~ 2 mag 
versus ~1 mag), may indicate that irradiation-induced en- 
hanced mass transfer is a factor at much shorter orbital peri- 
ods than for V344 Lyr. In addition, the amplitude of the ~0.5 
mag shoulder on the rise to superoutburst is roughly matched 
in our calculations; with enhanced mass transfer the shoulder 
amplitude is too great (> 1 mag). We consider it fortunate 
that we have matched the decay as well as we have since the 
scalings for the extrema in the limit cycle model, in particu- 
lar S m ax(«coid)» are most probably gross simplifications to a 
more complicated situation that is currently beyond our un- 
derstanding. 

8. CONCLUSION 

We have presented detailed time dependent calculations of 
the accretion disk limit cycle mechanism with application to 
the Kepler light curve of V344 Lyr. We find that the standard 
disk instability model is able to account for the mixture of nor- 
mal and superoutbursts with a minimal amount of fine-tuning. 
Therefore we argue that both the thermal-tidal model and the 
enhanced mass transfer model with their additional free pa- 
rameters, considered in detail by Schreiber et al (2004), are 
unnecessary. One can in some sense view the alternation of 
many successive short outbursts with a superoutburst as an 
extensive of the pattern of short and long outbursts seen in the 
longer orbital period systems such as U Gem and SS Cyg; for 
SU UMa parameters one has many more short outbursts be- 
tween two longs than in a system like SS Cyg, which exhibits 
LSLSLS... most of the time (Cannizzo & Mattei 1992). For 
our standard model parameters, we find that > 800 grid points 
are required for numerical stability, with a resulting pattern in 
which 5—7 short outbursts occur between two successive su- 
peroutbursts (for ft = 1/640). Having constrained a co id and 
cthot so as to produce the overall recurrence time scales of 
short outbursts and superoutbursts, there are several obvious 
deficiencies: the rise times of all outbursts are too fast, the 
superoutburst duration is too short, and the duration of the 
shoulder is too long and has an extra scalloped-out feature 
prior to the observed shoulder. However, the more physically 
motivated H98 scaling with T 0 = 12500 does a better job on 
making superoutbursts with minimal shoulders as observed, 
lasting d, and also longer viscous plateaus. Overall, our 
outbursts are too faint by ~1 mag for an intermediately in- 
clined disk. It may be that the assumed distance 620 pc for 
V344 Lyr (from Ak et al. 2008) is slightly too large. 
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